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ABSTRACT 


Problems of elastic stability of a space rocket under dynamic, non- 
conservative loads are presented. The spacecraft is treated as a slender, 
uniform cylindrical bar with free-free ends. Analysis of the frequency of 
the transverse vibration coupled with the longitudinal vibration under the action 
of time-varying thrust force based on Galer kin's method leads to Mathieu's 
equation. The elastic stability boundaries of the frequency ratio, longitudinal 
to transverse, versus the thrust are obtained. The dynamic buckling thrust 
of a uniform bar with various boundary conditions is presented. The aero- 
dynamic force induced from the oscillation of the rocket in supersonic speed 
on the critical thrust is investigated. A numerical illustration is given to 
show the relationship of the critical thrust versus flight speed at various 
flight altitudes. 
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DEFINITION OF SYMBOLS 


Symbol 

A 

A 

n 

a 

c 

c 

OO 

E 

I 

K 

L 

M 

m 

N 

P - P o 

p 

OO 

q 

R 

r 

n 

t 

t 

o 

U 


Definition 

Effective load carrying cross-sectional area 

= a \l2( 1 - cos 6 t )/ (0 t ) 2 ' 
no no 

Constant, defined by equation ( 11) 
slEA/m' 

Velocity of sound in air 

Modulus of elasticity 

Bending moment of inertia 

= P /(tt 2 EI/L 2 ) 
o 

Length of bar 
Mach number = U/c 

OO 

Mass per unit length of bar 
Axial force 

Thrust force function and its maximum value 
Atmospheric pressure 
Aerodynamic load per unit length 
Radius of bar 

= 9 /fii 

n 1 

Time variable 
Thrust buildup time 
Velocity of flight 


iv 



DEFINITION OF SYMBOLS (Concluded) 


Symbol 

u 


w 

x 

a 



Longitudinal displacement of bar 

nl A rrr\A fr* 0 £i tv* a ncuorco \ri f 1 An 

i.1 WX U1U.11UVW 111 VtV Vf* ww »*. «-»-« w ■w . — 

Transverse displacement 
Coordinate along bar 
- P o /mL»0/ K 

= jrRyP^M/mLfij 2 
= 1. 4, polytropy index of air 

= iTryRLp^/mLc^flj, aerodynamic damping factor 
Structural damping factor 
= n7rc/L, longitudinal frequency 

Parameter of Mathieu equation defined by equation (23) 

= t flj 
o 

Natural frequency of transverse vibration. 


v 
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Problems of elastic stability of a space rocket under dynamic, non- 
conservative loads are presented. The spacecraft is treated as a slender, 
uniform cylindrical bar with free-free ends. Analysis of the frequency of 
the transverse vibration coupled with the longitudinal vibration under the action 
of time-varying thrust force based on Galerkin’s method leads to Mathieu’s 
equation. The elastic stability boundaries of the frequency ratio , longitudinal 
to transverse, versus the thrust are obtained. The dynamic buckling thrust 
of a uniform bar with various boundary conditions is presented. The aero- 
dynamic force induced from the oscillation of the rocket in supersonic speed 
on the critical thrust is investigated. A numerical illustration is given to 
show' the relationship of the critical thrust versus flight speed at various 
flight altitudes. 


INTRODUCTION 


The means of finding the static buckling load of a slender column is the 
well known Euler's method. If the external forces have a potential, they are 
called conservative forces, in which case this method is applicable. If the forces 
are nonconservative, this method is not applicable [1], The theory of elastic 
stability of nonconservative problems is based on the investigation of the small 
oscillation of the system about its equilibrium position. Thorough discussion 
and many practical examples of this so-called dynamic method may be found in 
V. V. Bolotin's books listed as References 1 and 2. 

The thrust force and the aerodynamic load acting on a space vehicle have 
fixed directions relative to the vehicle and vibrate together with it when oscil- 
lation is taking place. Hence, these forces are nonconsei’vative. Determination 
of the elastic stability of the vehicle, while disregarding the local buckling and 
panel flutter, is the object of this study. The space vehicle may be treated 
simply as a uniform, slender cylindrical bar with both ends free. A thrust force 
is applied at the tail end along the body axis; this force is built up from zero to 
its maximum magnitude in a short time interval and then is kept constant. Two 
dynamic stability problems are presented here. 


0 


Dynamic Coupling of Longitudinal and Transverse Vibrations 

The velocity of the longitudinal stress waves traveling in a solid bar is 
equal to the speed of sound in the bar. In general, the frequency of the stress 
waves traveling back and forth in a bar is much greater than its natural frequency 
of transverse vibration; thus, the coupling is negligible. However, a rocket 
usually has great mass density per length and a small effective load-carrying 
cross-sectional area. Its longitudinal frequency may be of the same order of 
magnitude as its transverse vibration frequency. In this case the elastic coupling 
becomes significant. 


Elastic Stability of a Space Vehicle under the Combined Action 
of a Thrust and Aerodynamic Load 

For a space rocket of large cross-sectional area, the aerodynamic force 
and damping effect caused by its transverse oscillation in supersonic air flow 
may have some significance to its dynamic stability. Such information may be 
helpful to structural design engineers. 


ANALYSIS 


Coupling of Longitudinal and Transverse Vibrations 

The well known equations of transverse and longitudinal vibrations of a 
uniform bar and its boundary conditions of free-free ends are as follows [1, 3]: 


3*w 3 2 w 3 ( 3w\ . . 

+m w = ^ i N ^ / ) +q(x ’ 


( 1 ) 


N(x, t) = EA 

3x 


( 2 ) 


3 2 u 


3 2 u 


EA aV ' m e? 


(3) 


3 2 w 

3x 2 


x= 0 


3 2 w 

dx‘ 


x= L 


/ c) 3 w 

w 


x= 0 


S 3 w 

3x 3 


= 0 


7 x = L 


(4) 


2 



EA 


x=0 


(5) 



= -P(t) 



= 0 


The coupling of longitudinal vibration with the transverse vibration as given by 
the first term on the right-hand side of equation ( 1) can be obtained easily with 
the aid of Figure 1 in which N(x, t) denotes the tensile force along the axis of 
the bar. 


q(x,t) 


P ,,, f 4TITTm v 



X, u 


FIGURE 1. COORDINATE SYSTEM 


1* Solution of u(x, t) . - To eliminate the nonhomogeneous boundary 
condition of u given by equation ( 5) , we introduce a new variable u and write 
the solution of equation (3) in the form 


u(x, t) = u (x, t) - 


P(t) 

EA 


(x - x 2 /2L). 


( 6 ) 


From equations ( 5) and ( 6) , 



(7) 


Substituting from equation ( 6) into equation ( 3) yields 


^ . c 2 2% = m + _L (x _ x 2 /2Ll <£PM 

at 2 dx 2 mL EA ' /2L ’ dl? ‘ 


c 2 = EA/m . 


( 8 ) 
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Now, if we let 


u (x, t) = y u (t) cos , 
n= 0 n L 


which satisfies the boundary conditions given by equation (7) , the approximate 
solution of the unknown function u (t) may be determined using Galerkin's 
method. This involves substituting the assumed solution into equation (8) 
and then multiplying both sides of the equation by cos n7rx/L and integrating 
from x = 0 to x = L. This results in 


d 2 ^ 

^u 

n 

dt 2 


+ 9 2 u 
n n 


2L d 2 P( t) 
n 2 ^ dt 2 


n = 1, 2, .... 6 = mrc/L 
n 


(9) 


d 2 u 

~dF 


P(t) 

mL 


(10) 


The last equation represents the rigid body motion of the vehicle. 


Equation (9) indicates that the solution of u(t) depends on the second 
derivative of the thrust function P(t). It usually takes from one half to one 
second for a rocket engine to build up to its full power; then the thrust remains 
fairly constant until engine cutoff. The thrust buildup curve may be represented 
by a second order curve with t as the buildup time and P as its maximum 
thrust. Hence, let us assume 

P(t) = P [( 1 + ^ a) (t/t ) - ^ a(t/t ) 2 ] , 0 < a s 2, 0 < t < t ( 11 ) 

o 2 o 2 o o 


= P, t>t , 

o 

from which we obtain 


d 2 P(t) 

dt 2 


-P (a/t 2 ) rS(t) -S(t -t )], 
o o o 


( 12 ) 
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where S(t) is a unit step function. The thrust function and its second derivative 
is represented by Figure 2. 




FIGURE 2. THRUST FUNCTION AND ITS SECOND DERIVATIVE 


The solution of equation ( 9) may be obtained easily by using the Laplace transform , 
2LP A 

\ {t) = nVEA° oos ( V + V • 1 = V 1 13) 


where 


sin 0 t 


A = (a/T 2 ) n/ 2( 1 - cos r r )/r 4/ , 4> - tan 1 ~“7 , 

n o non n 1- cos 6 t ’ 


n o 


r 

n 


- ® n /^i» 


T = Sjt . 
O O 


In study of stability problems, the initial conditions can be disregarded; hence, 
we may drop out </> q from equation ( 13) and use P q to replace P(t) in dealing 
with the transverse vibration of the bar. Thus , we have 


u(x, t) = Yj 

n = 1 


2LP A 
o a 

n^EA 


„ x n?rx 

cos 6 t cos ~ 
n L 


- (x - x 2 /2L) P /EA 
o 



P /mDt 2 
o 

(14) 


and 


5 


(15) 


oo 2P A 

N(x, t) = - Y. — ~ ~ cos 6 t sin ^ 7 ^ - ( 1 - x/L) P . 

u , n 7 r n L o 

n= 1 


2. Solution of Transverse Vibration w(x,t) . -The solution of free 
vibration (N = 0 and q = 0) of a free-free beam is [4] 

00 

w(x, t) = Yj Cv -( x ) e^i* , 
i -1 1 1 

where Cb is a constant, j = \I - 1 , and v.(x) is the ith normalized vibration 


v.(x) = [Ch/3.x + cos /3.x - y.(Sh/3.x + sin /3.x) ] . 

1 \jh 1 iii i 

The symbols in the two equations above are defined as follows: /3. is the 
eigenvalue of the frequency equation 

cos j3 L Ch/3L = 1; 

S 2 is the natural circular frequency of free vibration 
fi . 2 = /3 . 4 El/m = (iS.L ) 4 ( 7 T 2 EI/L 2 )/ 7 r 2 mL 2 ; 

and y. is a constant 

1 

Ch/3.L - cos jS^L 
7 i = Sh/3.L - sin /3.L * 


The values of /3. and y. are 

/3.L = 4. 73000, 7. 8532, 10. 996, . . . , a £ ( 2i + 1) n 
y. = . 98250, 1. 00078, . 99997, . . . , a 1 . 


(16) 

mode, 

(17) 

(18) 
(19) 
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To solve equation (1) with q(x, t) = 0, we assume the solution in the form 


w(x, t) = Yj v.(x)f.(t) , (2°) 

i=l 1 1 


where f.(t) is an unknown function to be determined. After substituting w(x,t) 
from eq\iation (20) and N(x, t) from equation (15) into equation (1), multiplying 
the equation by v^v x) , and then integrating with respect to x from 0 to L, we 
obtain a system of infinite simultaneous equations: 


d 2 f, 


dt 2 


oo 

^ + n k 2(1 - Qa kk’ 2 2. v oos V )f i 


n= 1 kk 


ft, 2 y. (a +2 n cos 0 t) f. = 0, k = 1, 2, 
k i = i W n = i \i n 1 


i*k 


( 21 ) 


where 


a = P /mL 2 Sl/ 
o 


7T 2 


(W 


K = . 01972K 


K = P /(ti^EI/L 2 ) 

O 


a, . = ( ) 2 L f v [v.’ - (L - x)v!'] dx 
ki k J k l l 


( ft t /ft k ) 2 L [ Lv k ( 0) v.’ ( 0) + f (L - x) VjV.'dx ] 


p =-ffA (ftj/ft ) 2 L f v. ( ~ v" sin + v ' cos ) dx 
n^ n k k\mr l L i L 


n?rx 


y /- . 9 Jj r , , , n7rx . 

- aA (fty/ft, ) 2 — I v, 'v.' sin — — dx . 
n k mr -J k l L 
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The prime in the above equations denotes derivative with respect to x. Using 
the integration formulas given by Reference 5 we obtain 


■w-dg *i L 


. , 4P: % 3 

i + k, l k 


o 4 _ ^4 ~ Vi } + [1 " (_1) 1 (/?,/- /3, 4 ) l/" 7 ! 


k i 


( 22 ) 


k *i 


a.. = (fii/fijViS. L[i(y./3. L +6) - 4] . 
11 111 11 


(23) 


The approximate formula for fx may be written in the form 

ii 


« n s«(f.L|' [l-(-l) n iA ij 
ii 


4/3. 2 L 2 + n 2 7r 2 n 2 7r 2 


. 2/3. 2 L 2 + n 2 7T 2 

■i — J-T + 2 ‘ 


4/3. 4 L 4 + n 4 7T 4 


. (24) 


The above formulas give 

a n = 6. 151 a 12 = -9.211 

a 21 = . 4783 a 22 = 3 . 0 2 7 

Mi = hi = 7. 15 Ai = . 1423KA< . 
ii 

Because p decreases with 1/n 2 , let us consider only the first mode of both 
longitudinal and transverse vibrations. Then from equation (21) we obtain 

f ]2r 

-^y 1 + fij 2 ( 1 - aa^ - 2/Uj cos 0^) f^ = 0 . ( 25) 

Equation (25) may be rewritten in the form of the well known Mathieu equation, 
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where 


= ^i 2 ( 1 - aa u ) n = n/( 1 - or an) . 

Applying the stability boundary for the Mathieu equation given by Reference 2 
(pp. 24-28) , we have the critical frequency ratio, longitudinal to transverse: 


*1* = | n/ 1 - aa n ± n ' 

r l* = ^ n/i - Qa u +M n 2 / 3(1 - nan)' 

r i* = ~ n/ 1 - aay - 2^7 ( 1 - . 


(27a) 

(27b) 

(27c) 


The above equations for n = 1 and /q = . 1423KAj may be solved graphically, if we 
write these equations in the forms 


Aj = ± ( 1 - aa M - Irj 2 )/. 1423K 

(28a) 

A t 2 = ( 1 - aau - r t 2 ) ( 1 - oa u )/. 00675K 2 

(28b) 

Aj 2 = ( 1 - aau ~ r i 2 ) ( 1 " aa u )/ . 0405K 2 . 

(28c) 


Also, from equation (13) 

A t = a.\j2{ 1 - cos r i T 0 T/ ( r i T 0 ) 2 • 


(28) 


The procedure is as follows: 


( 1) Plot Aj vs rj from equations (28) and (28a), respectively, by using K 
as a parameter and note that aajj = . 1213K; a and t q are fixed, 

( 2) The intersections of these two curves give the upper and lower bounds 
of rj of the instability region for the given K. It is shown as the upper branch in 
Figure 3. 

(3) Plot A t 2 vs rj from equation (28b) and the square of equation (28) 
for given K,a, and t q . 

( 4) The intersection of these two curves gives the upper bound of the 
instability region of the lower branch shown in Figure 3. 

(5) Similarly, equations (28c) and (28) will yield the lower bound of 
the instability region of the lower branch. 

It is interesting to note that the thrust buildup time t Q (= t q/Dj) plays an 
important role in determining dynamic elastic stability of a space rocket. 


Dynamic Buckling Thrust and Frequency of Vibration 

Let us return to equation (25) and disregard the cos 0jt term in the 
equation; in other words, we will neglect the effect of the longitudinal stress 
we obtain 


+ (1 - aa n ) f t = 0 . (29) 

Therefore, the fundamental frequency of transverse vibration of a uniform beam 
under the influence of a thrust force is 


wave; thus. 



£2 = \Ti - aa^'Dj . 


(30) 


The critical dynamic thrust P * is the force that reduces the frequency of 
vibration to zero. It follows immediately that 


a* = i/a n or K * = P o */(7r 2 EI/L 2 ) = (/3 i L) A /ia n . 


(31) 
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For a free-free uniform beam, K* = 8. 23, and for beams with other end condi- 
tions, K* can be obtained easily with the aid of References 4 and 5. Formulas 
for various end conditions are given in Tables I and II. 

The Static Approach . -The critical thrust of a bar based on the static 
approach is actually the buckling inertia load caused by acceleration. It is assumed 
that the thrust force is increased gradually such that no transverse vibration is 
induced. In some cases, this approach fails to yield the critical inertia load, 
because under the given end conditions an equilibrium configuration of the beam can- 
not be attained without having transverse vibration. The critical values of thrust 
which causes inertia buckling of a clamped-free beam and a supported -supported 
beam, obtained from Reference 6, are .794 and 1.88 times ir^EI/L, 6 , respectively. 


Elastic Stability Under Aerodynamic Load 

It is a well known phenomenon that, as air flow passing around an elastic, 
slender body, energy is emitted from the body into the surrounding medium when 
the velocity of propagation of transverse wave of the body is greater than the 
velocity of the flow; and conversely, the energy is absorbed by the body from the 
medium. In the latter case the energy absorbed by the body has a tendency to 
increase its magnitude of vibration. This is commonly called self-excited 
motion or flutter. Our main concern here is not the panel flutter of the skin of 
the shell structure; we are interested in estimating how the high speed aero- 
dynamic flow affects the critical thrust. To this end, we shall avoid the compli- 
cated three-dimensional problem of an elastic shell in a potential flow such as 
that treated in Reference 1 (pp. 218 - 230). Our analysis will be made on the 
following assumptions: 

( 1) The deformation of the cross section of the hollow cylindrical bar is 
negligible; therefore, we may take the radial deformation of the rocket shell as 
the product of the transverse deflection of the free-free, uniform bar and cos 6 
with 0 measured from the plane of transverse motion around the cross section. 


w f (x, 6, t) = w(x,t) cos 0 . 


( 2) The bar is sufficiently long in comparison with its diameter so that 
the end effect may be ignored. 

{ 3) The linear approximation theory may be applied to determine the 
aerodynamic force exerted on the oscillating bar at supersonic speed. Let us 
denote U as the speed of the rocket or the velocity of air flow relative to the 
rocket, and by using the coordinates shown in Figure 4, we may write the 
linearized aerodynamic pressure in the form: 



TABLE I. CHARACTERISTIC FUNCTION AND EIGENVALUES OF VIBRATION OF A UNIFORM BEAM 
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TABLE n. CRITICAL THRUST OF A UNIFORM BEAM 
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( 32 ) 



Then the dynamic load per unit length, q(x, t) , in equation ( 1) is 


tt/ 2 fi-Ryp^ 

q = -2 f p cos 9 Rd 9 = 

-ir/2 C °° 


nv 


at 


9w 


- U 


ax j 


i). 


(33) 


N 


FIGURE 4. 


COORDINATE SYSTEM 


Stability Criteria. -Now, we proceed to solve equation ( 1) in which 
q(x, t) is substituted from equation (33) and N(x, t) is substituted from equation 
(15) by neglecting the longitudinal stres wave. Similarly, applying Galerkin's 
method, we obtain, by retaining only the first two terms, 




+ 2( ^a + y 121 


dfl 

dt 


+ fi 1 2 (l - aa u - /3b lt ) f t - +pbi 2 )l 2 = 0 


(34) 



+ 2(^ a 




dt 


+ fi 2 2 ( 1 ~ tta 22 " $^22) ^2 ~ ^2 2 ( aa 21 + ^ 2 l ) fl _ 0 , 


where 12 j and 12 2 are the first and second natural frequencies of free vibration, 
£ a and £ s are the aerodynamic and structural damping factors, respectively, 
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tfyRLPoo 

t = 

a 2mLc ’ 

oo 1 


and 


7rR-yp M 

OO 

mLfij 2 


(M = Mach No. = U/c ). (35) 

OO 


L 

b. . = (fii/fi, ) 2 L f v. v! dx = 4 [1 
ki 1 k ^ k l 
o 


( -i) i+k l(a 1 / V ^V(^-S l ‘). (3g) 


Note that the integration formula given by Reference 5 has been used to obtain 
the last expression, from which we calculate 


b n = 0 b 22 = 0 b^ = -9. 212 b 21 = 0. 1596 

Here, we point out that the rigid body modes have no effect on the elastic 
stability of a free-free bar. The proof is very simple. Let us add the rigid body 
motion f and rotation f to equation (20) to give 

-L\ U 


w(x, t)=f (t) + xf (t)+ y v. ( x) f . ( t) (20a) 

ri u . 1 1 

1-1 


Substituting the above equation, with 


N(x, t) = - P (1 - x/h) 
o 


and q(x, t) , given by equation (33) , into equation ( 1) and then multiplying the 
resultant equation by VjJx) , we integrate it from x = 0 to x = L. Using the 
integrals given by Reference 5, 


L 

f xv^(x) dx = 0 
o 


and 


L 

f v^(x) dx = 0 
o 


we note immediately that the equations obtained are identical with equation ( 34) . 

Substituting fj = Bje^*' and f 2 = B 2 e^ into equation ( 34) and equating the 
determinant of the coefficients of Bj and B 2 to zero yields the characteristic 
equation 
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where 


a = 1 
o 

a 1 -2efi 1 e = 2(4 + 4 ) 

as 

a 2 = aj 2 [(i + r + e 2 ) - o(a H + r a 22 ) ] r = ( S2 2 /S2 A ) 2 = (/3 2 //J t ) 4 = 7. 5985 
a 3 = a 4 3 [i +r - o^a^ + r a 22 )]e 

a 4 = [1 - cHau + a 22 ) + a 2 (aua 22 - aj 2 a 2 j) - (a^b^ + a 2 jb y)otP 

- bcb^/S 2 ] . 


We recall that the Routh-Hurwitz stability criteria are as follows: 


(a) aj > 0 ; or e > 0 


(38) 


(b) 


a l a 3 


a a 2 

o ^ 


n * i^+l+r 

> 0 ; or a < 2 — 

a ll + r a 22 


“ .294 


(39) 


(c) 


a l a 3 0 

a Q a 2 a 4 
0 aj a 3 


> 0 or Ao 2 - Bo: + C > 0, 


(40) 


where 


A = (a u + ra 22 )/4r - (a 11 a 22 - aj 2 a 21 ) = 4. 1512 



( 1 +r + e 2 ) (a 11 + ra 22 ) - [an + a 22 - ( a 12 b 21 + a 21 b 12 )/3 ] 


= 1. 9263( 8. 5985 + e 2 ) - ( 9. 194 - 5. 8763/3) 


C = ( 1 + r) ( 1 + r + 2e 2 ) /4r - ( 1 - b 12 b 21 /3 2 ) = . 2829(8. 5985 + <r 2 ) 


- ( 1 + 14. 7 0 2 8/3 2 ) 
(d) d 4 > 0 . 


Condition (a) is satisfied if the damping factor is positive; and (d) is satisfied 
for all values of a; (b) requires that a* be smaller than . 294 or K* < 14. 9; 
and (e) can be used to determine the critical thrust for a given ( 3 . 

Frequency of Transverse Vibration . -Let us set the damping factor equal 
to zero (aj = a 3 = 0) in equation (37) from which we may readily solve for the 
frequency of vibration by using the thrust P or the aerodynamic factor j3 as a 
parameter. A one-term approximation without taking into account aerodynamic 
effect is given by equation ( 30) . Note that the coalescence point of these curves, 
by definition, is the limit of stability of motion. 


NUMERICAL EXAMPLE 

Let us consider a large rocket with the following structural data: 

R = 4. 572 m (15 ft) , 

L - 106. 68 m (350 ft), 

Weight = mgL = 2. 72 x 10 6 kg ( 6 x 10 6 lb) , 

El = 1. 4632 x lO 10 m 2 -kg ( 50 x 10 12 in 2 -lb) , 
from which we calculate 

fij = 4. 66 rad/sec and tt 2 EI/L 2 - 12. 7 x 10 6 kg (28 x 10 6 lb). 

The values of t, and 1 3 given by equation (35) are given in Table III. 

3 , 

By using equation (40) the critical thrust versus the flight speed of the 
rocket is plotted in Figure 5 with flight altitudes at 3048 m ( 10, 000 ft) , 6096 m 
(20,000 ft) , and 9164 m (30,000 ft) , respectively. These curves indicate that 
structural damping has little effect on the critical thrust. 
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TABLE HI. VALUES OF £, AND p 

a 


Altitude 

‘a 

P 

3048 m (10,000 ft) 

.0192 

. 0236 M 

6096 m (20,000 ft) 

.0137 

. 0158 M 

9164 m ( 30, 000 ft) 

.0090 

.0102 M 

15,240 m (50,000 ft) 

.0036 

. 00394 M 

30,480 m (100,000 ft) 

.00033 

. 00036 M 

42,672 m (140,000 ft) 

.000034 

. 000068 M 


The variation of vibration frequency, which is represented by equation 
( 37) with aj = a 3 = 0, is illustrated by Figures 6 and 7. In Figure 6 the frequency 
ratio, frequency with thrust to frequency without thrust, is plotted against P q 
with p as a parameter, while in Figure 7 it is plotted against ft with P Q as a 
parameter. Note that the one-term approximation given by equation ( 30) is also 
shown in Figure 6. For a given rocket the value of ft is proportional to Mach 
number at a constant altitude. This is also shown in the figures. 


CONCLUSIONS 


( 1) It has been proven that the rigid body modes do not affect the 
elastic stability of a uniform bar with free-free ends. 

(2) From a structural design point of view, the critical thrust is more 
than four times greater when the rocket is treated as free-free with an axial 
thrust than when the rocket is treated as supported-supported with a vertical, 
conservative force. 

(3) Figure 3 indicates that the region of elastic instability decreases 
with the thrust buildup time of the rocket. 

( 4) The critical thrust decreases rapidly at high supersonic speed and 
low flight altitude. 
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( 5) Structural damping does not have significant effect on the critical 

thrust. 

(6) One must keep in mind that the aerodynamic load given by equation 
(33) is a very crude approximation, but that it leads to the simplest formulation 
to serve the purpose of our study. Certainly, one can find many other approaches 
based on more realistic assumptions. Discussion of these approaches is beyond 
the scope of this report. 
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Mach Number 


FIGURE 5. CRITICAL THRUST VS FLIGHT SPEED 
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FIGURE 6. VARIATION OF FREQUENCY VS THRUST 
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